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Abstract. Systems with absorbing (trapped) states may exhibit a nonequilibrium phase transition from 
a noise- free inactive phase into an ever-lasting active phase. We briefly review the absorbing critical phe- 
nomena and universality classes, and discuss over the controversial issues on the pair contact process with 
diffusion (PCPD). Two different approaches are proposed to clarify its universality issue, which unveil 
strong evidences that the PCPD belongs to a new universality class other than the directed percolation 
class. 
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1 Introduction 

The absorbing phase transition (APT) which may occur in 
a system with absorbing states has been an active research 
topic in nonequilibrium statistical mechanics with possi- 
ble applications to a wide range of areas in physics, chem- 
istry, biology, and sociology |l|2|3|4j . Microscopic states 
into which probability is collected are called absorbing. 
Once a system is trapped into absorbing states, it can- 
not escape. Absorbing states sometimes form an absorb- 
ing space of many states around which a system wanders 
forever, but does not escape out of the absorbing space. 
Such an example is the state with only one particle dif- 
fusing in physical space for the pair contact process with 
diffusion (PCPD); see Sect. [3 

In the study of the APT, two different quantities may 
serve as order parameter^. One is the steady state den- 
sity of the activity (outside of the absorbing states) in the 
thermodynamic limi10, say p whose definition varies from 
system to system, and the other is the survival probabil- 
ity, say P s , with which a system docs not get trapped into 
absorbing states forever. Since a finite number of activity 



1 Although the APT has nothing to do with the order- 
disorder transition, the jargon "order parameter" has been 
used to name the indicator of the phase transition which takes 
zero in one phase and non zero in the other. 

2 For systems with finite volume of the configurational state 
space, the presence of the absorbing state always lets a system 
fall into it eventually [5]. So, as in the equilibrium statistical 
mechanics, the thermodynamic limit (infinite extension of the 
configurational state space) is indispensable to study a non- 
trivial APT. 



(or local activity) means zero activity density in the ther- 
modynamic limit but non zero survival probability, p = 
in general does not imply P s — although the reverse is 
always true. In a sense, p (P s ) is the order parameter for 
the "macroscopic" ( "microscopic"0) APT. Hence, these 
two quantities might pinpoint different transition points 
in principle, see Sect. [3] for an example. However, with 
proper initial conditions and definitions for each order pa- 
rameter, these transitions points may coincide in general. 

Besides the difference mentioned above, nature of the 
phase transition described by these order parameters can 
be different even if they locate the same critical point. 
To clarify this point, consider the branching process with 
spontaneous death (A — > 2A with rate a and A — » 
with rate A). The absorbing state of this model is a state 
without A. This problem is the linear one step process in 
Rcf . [5 , which can be solved easily and one may find the 
survival probability as 



f if A > cr 
P '= ^ ifA<„ 



(1) 



The term "microscopic" is employed because the survival 
probability has been mainly analyzed in practice when the ac- 
tivity is initially localized in an infinite lattice. In this case, 
the macroscopic density p remains zero for all time by defini- 
tion though the number of activity may increase indefinitely 
with finite probability. Hence, the meaningful activity density 
p in the thermodynamic limit should be examined with a finite 
initial density. 
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On the other hand, the mean number of particles (activ- 
ity), say (n), behaves as 

^L = {a-\)(n)^(n)= m e^\ (2) 

with the initial value no- Both quantities are singular at 
A = a, but P s increases continuously unlike (n) which 
shows a sharp jump at the transition point. Hence in the 
study of the APT, P s and p should be studied indepen- 
dently and both quantities are important in understanding 
the APT. 

In what follows, however, we restrict ourselves to the 
study of the "macroscopic" APT, or the scaling behav- 
ior of p and quantities related to it like the correlation 
functions. Like equilibrium phase transitions, the APT is 
characterized and classified by the critical exponents; the 
order parameter exponent f3, correlation length exponent 
UX; relaxation time exponent U\\, and so on, which are 
defined as 

P~(Pc-pf, ^\P-Pc\- U \ T~\p- Pc \- V \ (3) 

where p is the external tuning parameter with p c to be 
the critical point, £ is the correlation length, and r is the 
relaxation time. One can define other critical exponents 
but with the aid of the scaling ansatz most of critical ex- 
ponents can be deduced from the above exponents [6]. 
For example, the density decays as t~ s at criticality with 

To categorize universality classes according to the crit- 
ical exponents and to understand what properties combine 
different systems into the same universality class are the 
main goals in this field. Some understanding has emerged 
from the numerical and analytical studies. Section[5]briefly 
summarizes the well-established universality classes, such 
as the directed percolation (DP), the directed Ising (DI), 
the parity conserving (PC) classes and so forth. The last 
10 years have witnessed intensive discussion and hot de- 
bate on a simple but very elusive interacting particle sys- 
tem, the PCPD. In Sect. [31 we will critically review on the 
issue regarding its universality class. To settle the contro- 
versy, two different approaches have been proposed by the 
authors, which is the subject of Sect. 2] We draw a con- 
clusion in Sect. [5] 



2 Universality classes 

The simplest non-trivial model which shows an APT might 
be the contact process (CP) which is an interacting hard 
core particle system on a d dimensional lattice with the 
creation of a particle by a neighboring particle and the 
spontaneous annihilation (A — > 2 A, A ~ ► 0) [7]. The par- 
ticle vacuum state is the only absorbing state in which 
the system cannot escape by the prescribed rules. Hence 
the order parameter is the density of occupied sites (or 
particles). The CP shares the critical behavior with the 
directed percolation (DP) the preferred direction of which 
is interpreted as the time direction of the CP and the open 



channel of which as a particle 2 . The DP has the rapidity- 
reversal (or time-reversal) symmetry [618] which associates 
the microscopic APT with the macroscopic one, that is, 
which renders P s to scale equivalently to p. A nice illustra- 
tion of the connection in the context of the bond directed 
percolation can be found in Ref. [2]; see also Ref. [4]. 

After extensive numerical studies regarding the uni- 
versality class of the APT, it has been conjectured^ that 
the APT model with a single absorbing state should be- 
long to the DP class if symmetry or conservation is not 
involved [819110] . The robustness of the DP class extends 
to the systems with infinitely many absorbing states like 
the pair contact process (PCP) [11], at least in its sta- 
tionary property. In the PCP, the creation and annihila- 
tion of particles are only mediated by two particles which 
form a nearest neighbor pair on a lattice (2A — > 3A, 
2 A — > 0). So any configuration devoid of a pair is an ab- 
sorbing state and the number of absorbing states increases 
exponentially with system size. Clearly, the order param- 
eter should be the pair density not the particle density 
(auxiliary field density) which is nonzero at stationarity 
irrespective of the phase. The robustness of the DP class 
suggests that the absence of symmetry or conservation of 
the order parameter should render the system to belong to 
the DP class, irrespective of the existence of the auxiliary 
field associated with infinitely many absorbing states. Al- 
though the stationary property of the PCP conforms with 
the DP scaling, there is still vivid discussion regarding 
the dynamic scaling of its spreading |12ll3ll4ll5j . which 
is beyond the scope of this paper. 

Other universality classes have been found by adding 
symmetry or conservation in dynamics. The directed Ising 
(DI) class involves the Z2 (Ising) symmetry in dynam- 
ics, the evolution operator of which is invariant under 
the Z2 symmetry operation [16117118] . Naturally, the typ- 
ical DI systems include two equivalent absorbing states. 
The DI scaling also applies to systems with two equiv- 
alent groups of multiple absorbing states [19120121] , The 
conservation of the particle number of modulo 2 reveals 
another universality class called as the parity conserving 
(PC) class [22123] . which coincides with the DI class in 
one dimension. In higher dimensions, both classes can be 
described by the trivial mean-field theory. Recently an- 
other universality class (generalized voter class) has been 
examined, which also coincide with the DI class in one 
dimension [53]. It is worthy to note that the DI (and PC) 
class returns to the DP class immediately with the intro- 
duction of a symmetry breaking field or a conservation 
breaking dynamics [25126127128] , 

There had been an attempt to find a new universality 
class by studying models with higher symmetry than Z2 
or mod(g) conservation with q > 2. However, all models 
so far show a trivial critical behavior even in one dimen- 
sion, in that the absorbing phase is found to be always 
unstable against the dynamics increasing activity of the 
order parameter. 

All the systems explained up to now can be described 
by the single component order parameter. Richer behavior 

4 This statement is termed as the "DP conjecture" . 
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is anticipated when multi species are involved. For exam- 
ple, the interaction of the order parameter with a con- 
served field triggers a different type of universal behavior 
depending on the activity of the conserved field |29|30j . 
In fact, any system involving multi-particle reactions can 
be interpreted as a multi-species particle system. One can 
map the PCP to the multi-species model by identifying a 
pair as a particle of one kind and a single isolated particle 
as a particle of another kind. More general cases will be 
discussed in subsequent sections. 



3 Pair contact process with diffusion 

The pair contact process with diffusion (PCPD) is an ex- 
tended model of the PCP with hopping of particles allowed 
for. To be specific, the dynamics of the PCPD in d dimen- 
sions consists of hopping, pair annihilation, and creation 
by a pair, which is symbolically summarized as 



AAQ\ 
UAj 



AA -> 00 with rate 
> AAA with rate 



P 

d' 
l-p 



2d 



(4a) 
(4b) 
(4c) 



where A (0) stands for an occupied (a vacant) site on a 
lattice and < p < 1. The PCP corresponds to the above 
rules with D = which make any configuration without 
two particles in a row (a pair) absorbing. Due to the dif- 
fusion, however, a state without a pair but many isolated 
particles is not absorbing any longer. Only both the par- 
ticle vacuum and the state with only one particle in the 
whole system are absorbing. Since the particle density of 
both absorbing states of the PCPD is zero in the thermo- 
dynamic limit, it can play the role of the order parameter 
in contrast to the PCP case. Needless to say, the pair den- 
sity may also serve as an order parameter. 

Due to the lack of a process to eliminate a single par- 
ticle (no single particle reaction) without particle colli- 
sions, the conventional survival probability and the den- 
sity might locate different transition points. To elucidate, 
consider the PCPD in higher dimensions than 2. If initially 
two particles are located somewhere in an infinite lattice 
(just outside of the absorbing space), the survival proba- 
bility is always finite because of the nonrecurrence of the 
random walk even for the case of p = 1. That is, the sur- 
vival probability predicts the absence of the "microscopic" 
absorbing phase. However, the macroscopic critical point 
should be located at finite p. The reason is as follows: The 
mean density p of particles in d dimensions satisfies the 
(exact) equation, 



^ = (1 - 3p)p p - (l-p)p t , 



(5) 



where p p (p t ) means the pair (triplet) density. lip > i, the 
steady state value of the pair and triplet density should 
be zero, which is clear by Eq. ((5|). If p were not zero in 



the steady state, macroscopic number of pairs should be 
formed by the diffusion, which is contradictory to the ob- 
servation that p p — 0. Hence p should approach to zero 
if p > | and the critical point should be not larger than 
^ and the "macroscopic" absorbing phase is present in 
any dimension^]- In the above argument, we assume the 
existence of the steady state even in the thermodynamic 
limit. 

As the above consideration reveals, the single-particle 
diffusion plays a crucial role in changing the nature of the 
conventional microscopic APT, which has nothing to do 
with the DP for higher dimensions than 2. However, this 
does not resolve the controversy regarding the universality 
class of the PCPD. First, the difference of the microscopic 
APT does not guarantee that of the macroscopic one. A 
good example is the PCP whose macroscopic APT is char- 
acterized by the DP scaling but whose microscopic APT is 
known to be nonuniversal |12|13|14|15j . Second, the main 
issue is not any dimensional PCPD but one dimensional 
PCPD where it is not fully clear that the microscopic APT 
is equivalent to the macroscopic APT |31|32j . 

In fact, the difference between the PCPD and the DP 
is well appreciated for two or more dimensional systems. 
First consider the mean field equation. The mean field 
equation for the PCPD can be found by replacing p p and 
pt with p 2 and p 3 , respectively in Eq. |5j) which reads 



^ = (1 

dt y 



3 P )p 2 -(l- P )p 3 



(6) 



The mean field critical exponents are (3 — 1 and 6 = 
f3/v\\ = | [331 34] . which are different from those of DP 
(J3 = vn =1) [2J. Since the upper critical dimension of the 
PCPD is believed to be 2 [35], for most physically relevant 
cases (d > 2) the PCPD does not belong to the DP class. 

For the one dimensional PCPD, however, the numeri- 
cally estimated critical exponents of the PCPD are so sim- 
ilar to those of the DP that the possibility for the PCPD 
to belong to the DP class has been raised [36l37j . Interest- 
ingly, the critical exponent^ 6 which describes the density 
decay with time at criticality has floated from ~ 0.28 [38] 
to less than 0.185 |37 with time, which is due to the strong 
corrections to scaling. For comparison, the numerical value 
of S B p is ~ 0.15946 [33]. 

An argument in favor of the DP scenario was suggested 
by Hinrichsen [37j . The starting point of the argument is 
the numerical observation that the dynamic exponent z 
is smaller than 2 which is the dynamic exponent of the 
random walk. Hence, if coarse graining is performed ac- 
cording to the PCPD dynamic exponent, the diffusing iso- 
lated particles will stop moving asymptotically and the 



5 To resolve the difference between the microscopic and 
macroscopic APT, a new definition of survival of the system 
has been suggested: The system without a pair is considered 
as (temporarily) absorbing. With this definition, it has been 
shown that two order parameters exhibit the APT's at the 
same transition point at least in one dimension [31132] . 

6 From now on, by critical indices without subscript are al- 
ways meant those of the PCPD. 
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long time behavior of the PCPD should be same as that 
of the PCP which belongs to the DP. 

We would like to make some comments as to this ar- 
gument. To begin, it is not at all clear why the dynamics 
of the wandering isolated particles is decoupled asymptot- 
ically from that of the active clusters. If not, the dynamic 
scaling of the isolated particles should be affected by the 
complex environmental geometry of active clusters. There- 
fore there is no ground for the belief that the diffusion of 
isolated particles remains governed by the random walk 
dynamic exponent zrw = 2. If the system scales in one 
way as a whole (not decoupled), the coarse graining argu- 
ment does not lead to the zero diffusion constant of the 
isolated particles. Second, the Hinrichsen's argument set 
the upper bound for the dynamic exponent z < 2 in or- 
der to be self-consistent. Based on this, one can comment 
on the two dimensional PCPD which seems definitely not 
in the DP class. This leaves us only one option to take 
z = 2 for the two (and higher) dimensional PCPD even 
without any logarithmic correction. However, the upper 
critical dimension of the PCPD is believed to be 2, where 
the logarithmic correction is expected and the numerical 
results seem to support its presence [35]. Therefore, even 
if one may accept the decoupling of two different fields in 
the PCPD, his argument seems not working at two dimen- 
sions. 

Regardless of the universality issue, one can ask why 
the PCPD has such strong corrections to scaling which 
are the main obstacle in numerical study. The long term 
memory effect was suggested as a possible origin of the 
strong corrections to scaling though it is not clearly an- 
swered how the finite mean life time from the life time 
distribution P(t) ~ r~ 2 - 2 can trigger such strong correc- 
tions to scaling or even a new scaling deviated from the 
DP |31| . Only a possible scenario was contemplated that 
the interacting theory may force the system to flow into a 
non-DP fixed point even with finite mean life time, reflect- 
ing on the similar result found in the Levy flight DP sys- 
tems |4QI41j . Due to the fact that the memoryless case of 
the generalized PCPD clearly shows the DP scaling [3"T] , 
it is certain that the memory effect plays an important 
role in determining the universality class of the PCPD. 
Another reason of believing the role of the memory comes 
from the model tpl2 (for the definition, see below) and 
similar varieties of models with rules nA — > (n + m)A 
and I A — > (I — k)A with n > I and m, k > which are 
numerically shown to belong to the DP [42] , Same as the 
PCPD, the tpl2 (n = 3, 1 = 2, m = 1, k = 2) with zero dif- 
fusion constant has infinitely many absorbing states and 
belongs to the DP class [43] . In contrast to the PCPD 
(n = 2, 1 = 2, m = 1, k = 2), the diffusion turns out not 
to affect the universality class of the tpl2. In Ref. [43] , 
it is argued that the effective lack of the memory in the 
tpl2, which is clearly seen from the space-time configura- 
tion Fig. [H renders any model with n > I to belong to the 
DP class. 

Another complicated feature of the PCPD has arisen 
in the analytical study. The formulation of the field the- 
ory based on the single field turned out impossible and it 



Fig. 1. (Color online) Space-time configuration of the tpl2 at 
criticality with a sparse initial condition. The time flows from 
top to bottom and the cross section by a horizontal line is the 
space configuration at specified time. The periodic boundary 
conditions are employed. Effective lack of the memory (see the 
text) means the seemingly absence of the branching process 
3A — > 4A due to the low probability of forming a triple before 
a pair is annihilated. 

is argued that the nonperturbative treatment still cannot 
cure this failure [H]. The lesson from the study is that 
at least two independent field should be included in the 
proper action to understand the PCPD in the field theo- 
retical framework. The proper field theory which is local 
in space and tim^Zl, if exists, should take two independent 
fields into account. In principle, one can write down the 
Langevin equation equivalent to the PCPD using two in- 
dependent fields [3S], but the proper analytical tool does 
not seem to be at hand. 



4 How to tackle the problem 

The one dimensional PCPD seems very difficult, if not 
impossible, to tackle directly. In the numerical front, the 
strong corrections to scaling prohibits researchers from 
measuring the critical exponents accurately. In the an- 
alytical front, the failure of the field theory by the single 
field requires an ingenious treatment of the problem. 

In this section, we will try to convince readers that the 
PCPD does not belong to the DP class. First, introducing 

7 The effective field theory might be described by a sin- 
gle field with nonlocal interaction, which in principle can be 
achieved by integrating out some of the fields in a local action. 

8 The most recent numerical studies still do not provide a 
conclusive evidence [46147] . 
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Fig. 2. (Color online) Semi-logarithmic plot of the pair density 
multiplied by t s ° p with fop = 0.1594 vs time for the PCP 
with biased branching. In the active (absorbing) phase, the 
curve veers up (down) and at criticality the flat straight line is 
expected, which renders the estimation of the critical point as 
p c = 0.077 168(1) with error in the last digit. 

the bias, we will argue that the failure of the field theory 
by the single field is also generic even in one dimension. 
Second, studying the crossover model from the PCPD to 
the DP by introducing the single particle dynamics, the 
difference of these two classes will be clarified. 



4.1 Effect of Biased Diffusion 

As discussed at the end of Sect. El the PCPD is supposed 
to be described by two independent "elementary excita- 
tions" , that is, the isolated particle-field and the pair-field. 

Even if we find the proper field theory for the PCPD, 
it might not resolve the controversy on the one dimen- 
sional PCPD. For example, what if the fixed point of the 
PCPD turns out not to be reachable by the perturbation 
expansion to the one dimensional system just as that of 
the branching annihilating walk (BAW) with even num- 
ber of offspring [48] ? However, before worrying about the 
above scenario, we must take a first step to answer the 
most elementary question whether it is absolutely neces- 
sary to employ two independent fields in the one dimen- 
sional PCPD, in contrast to the Reggeon field theory of 
the DP which needs a single field. The answer may not 
resolve the controversy definitely, but should be regarded 
as a big step in understanding the difference between the 
PCPD and the DP class. 

Recently, the present authors have suggested to check 
if the one dimensional PCPD should have two independent 
relevant fields via adding a bias in diffusion [49 34J . We got 
some hints from the study of the two species annihilation 
model A + B — > which shows different decaying behavior 
in the presence of the relative bias between A and B from 
that without bias [50 51]. 

Since the DP has only one relevant field, the bias should 
be gauged away by the Galilean transformation. Besides, 
lots of complicated models which belong to the DP class 



also turn out to be robust against the bias |49|43j . Im- 
plementing the relative bias between relevant fields and 
observing how the system responds can be a litmus test 
to check numerically whether the single field is enough to 
describe the system. In this context, we introduced the 
driven pair contact process with diffusion (DPCPD) and 
studied it [49|34j . 

Before delving into the DPCPD, we should emphasize 
that the presence of two "independent" fields is not suf- 
ficient to observe the change of critical behavior by the 
bias. Although the PCP is described by two independent 
fields, we expect that the branching bias which triggers 
the relative bias between pairs and isolated particles does 
not alter the critical behavior because the PCP belongs to 
the DP class. In Fig. [21 we showed how the pair density 
behaves around the critical point of the PCP with biased 
branching. As in the PCP, p is the pair annihilation prob- 
ability and 1 — p is the creation probability. The branching 
bias is given such that a created particle is always placed 
at the right neighbor of the pair once that site is vacant. 
Figure [2] locates the critical point p c as p c = 0.077 168(1) 
by exploiting the power law behavior of the pair density, 
say p p , at criticality. As can be seen, the critical decay ex- 
ponent still takes the DP value (<5dp), which implies that 
the relative bias is irrelevant in the PCP. 

Now let us continue the discussion on the DPCPD. In 
Ref. [49], the DPCPD shows clear distinction from both 
the DP and the PCPD scaling in one dimension, which 
was confirmed again by the cluster approximation along 
with the coherent anomaly method analysis [34] . This re- 
sult is not expected if there is only one relevant field in 
the PCPD. Hence the field theory of the PCPD should 
be different from the Reggeon field theory even in one di- 
mension, which provides a strong evidence of a non-DP 
scaling in the PCPD. 

Another interesting feature of the DPCPD is the di- 
mensional reduction. The one dimensional DPCPD is nu- 
merically found to have the same critical behavior as the 
two dimensional PCPD [49]. The dimensional reduction 
by the biased diffusion is actually reported in different ar- 
eas in nonequilibrium statistical physics. One example is 
the two species annihilation model, A + B — > 0. When 
there is (no) relative bias between A and B, the density 
decays as t"( d+1 )/ 4 (t~ d / 4 ) wh en d < 3 (d < 4) and as 
< _1 when d > 3 (d > 4) [50 5 lj . Another example can be 
found in the study of the self organized criticality. It is 
rigorously proved that the upper critical dimension of the 
directed sand pile model is 3 unlike its undirected version 
whose upper critical dimension is 4 52J. 

4.2 Learning from crossover scaling 

Another evidence that the PCPD does not belong to the 
DP was provided by the study of the crossover from the 
PCPD to the DP [53] . Along with the dynamics in Eqs. §4$ , 
we introduce the single particle annihilation/creation 

A^0, Z}^^AA, (7) 
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where < q < 1. li w ^ 0, the system shows the DP 
scaling behavior [53] , 

What can we expect about the scaling behavior near 
the PCPD critical point for finite w if PCPD does be- 
long to the DP class? From the rigorous study about the 
crossover "0 model of the BAW with one offspring [43] 
based on the stochastic equivalence shown in [S3], the 
"crossover" between the same universality classes is ex- 
pected to have two characteristics. First, the phase bound- 
ary is expected to meet the PCPD point linearly. Sec- 
ond, the critical amplitudes which is defined as p(t)t 5op 
for cases with finite but small w should collapse with the 
critical amplitude of the PCPD at w = 0. Our recent nu- 
merical results show a nontrivial crossover exponent from 
the PCPD to the DP and the phase boundary approaches 
the PCPD point in a singular (nonlinear) way [53]. This 
provides another strong evidence that the PCPD is differ- 
ent from the DP. 

In fact, one should be cautioned in interpreting the 
singular behavior of the phase boundary. The linearity of 
the phase boundary may become complicated by the non- 
trivial singularity arising between the PCP and the DP 
crossover [43], which corresponds to the dynamics mod- 
eled by Eqs. (j4j) and ([7]) with D — 0. However, this sin- 
gularity has nothing to do with the universality class. In 
a sense, the model with w = 0, i. e., the PCP, is patho- 
logical in that the configuration volume occupied by the 
absorbing states is macroscopic, which is not the case for 
finite m °l Hence there is an inherent singularity, actually 
discontinuity of the particle (auxiliary field) density, close 
to w = 0, which is reflected by the nontrivial crossover ex- 
ponent although the PCP belongs to the DP class [43]. If 
we introduce 3 A —> rather than Eq. {7J, we reproduced 
two characteristics of the "crossover" between models in 
the same universality class; see Fig. [3] The results sum- 
marized in Fig. [3] are rather easily conceivable because the 
operator corresponding to 3A — > is irrelevant in the RG 
sense and moreover there is no singularity of the auxiliary 
field density near w — 0. In this context, the single particle 
branching/annihilation introduced to the PCP is relevant 
because this operator changes the structure of the absorb- 
ing configurations, which is manifest by the singularity at 
w = [43], though it does not change the universality 
class. 

Since the volume of the absorbing states in the config- 
uration space for the PCPD is zero in the thermodynamic 
limit, there is no singularity arising from the auxiliary 
field near the crossover from the PCPD to the DP. Hence 
the reason of the singularity in the PCPD-DP crossover 
should be understood in another context. The easiest an- 
swer may be that the PCPD does not belong to the DP. 
Of course, there could be a hidden unknown reason for 
the PCPD to lead to a singular crossover behavior and its 
possibility cannot be fully excluded. Hence we studied a 



We use the quotation mark to mention the "crossover" be- 
tween models in the same class because the "crossover" is not 
truly a crossover in convention. 

10 When w is finite, the only absorbing state is the particle 
vacuum. 
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t 

Fig. 3. (Color online) Log-log plots of the particle density 
(upper curves) and the pair density (lower curves) at criti- 
cality vs t for the PCP "crossover" model with the dynam- 
ics 3A — > which occurs with probability w. The critical 
points are 0.077 0905(5), 0.07651(1), and 0.07594(1) for w = 0, 
w = 5 x 10 -4 , and w — 10 -3 , respectively. The difference 
among cases with w = and finite w is hardly seen. Inset: The 
phase boundary near w = 0. The line meets the vertical axis 
linearly. 

similar model to the PCPD whose non diffusing counter 
part has infinitely many absorbing states and which is 
known to belong to the DP class. It is found that the 
crossover model based on the tpl2 have two properties 
of the "crossover" among models belonging to the same 
universality class [J3]. We are unable to conceive of a 
mechanism that would explain these observations, while 
maintaining the PCPD and tpl2 in the same universality 
class; a more natural interpretation is that the two models 
belong to different classes, so that the PCPD does not fall 
in the DP class. All together, our study of the crossover 
scaling strongly suggests that the PCPD does not belong 
to the DP class. 



5 Conclusion 

Up to now, we discussed about the hotly debated issue of 
the pair contact process with diffusion (PCPD). Among 
many scenarios proposed at the first stage of the contro- 
versy [33], only two seem to have survived; does the PCPD 
belong to the directed percolation class or form a differ- 
ent universality class from any other known one? In this 
paper, we gave evidences in favor of the second scenario. 

At first, the fact that the biased diffusion drastically 
changes the critical behavior of the PCPD suggests that 
the one dimensional PCPD should be described by the 
two independent relevant "elementary excitations" which 
is not the case of the Reggeon field theory (DP). Second, 
the nontrivial crossover scaling arising between the PCPD 
and the DP strongly suggests that the PCPD should not 
belong to the DP class. 
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Although the main interest of this paper is the PCPD, 
the methods employed to clarify the universality issue 
of the PCPD (the role of the biased diffusion and the 
crossover) are generally applicable to many other systems, 
some of which are currently under investigation. 
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